EPJ manuscript No. 

(will be inserted by the editor) 



Static observables of relativistic three-fermion systems with 
instantaneous interactions 



o 
o 



in 

(N 
O 

o 

a.: 



HISKP-TH-06/02 

Christian Haupt. Bernard Metsch and Herbert-R. Petry 

Helmholtz-Institut fiir Strahlen- und Kernphysik, Nufiallee 14-16, D-53115 Bonn, Germany, e-mail: haupt@itkp.uni-bonn.de 
February 16, 2006 

Abstract. We show that static properties like the charge radius and the magnetic moment of relativistic 
three-fermion bound states with instantaneous interactions can be formulated as expectation values with 
respect to intrinsically defined wavefunctions. The resulting operators can be given a natural physical 
interpretation in accordance with relativistic covariance. We also indicate how the formalism may be gen- 
eralized to arbitrary moments. The method is applied to the computation of static baryon properties with 
numerical results for the nucleon charge radii and the baryon octet magnetic moments. In addition we make 
predictions for the magnetic moments of some selected nucleon resonances and discuss the decomposition 
of the nucleon magnetic moments in contributions of spin and angular momentum, as well as the evolution 
of these contributions with decreasing quark mass. 
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1 Introduction 

Static observables of bound state systems in field theo- 
retic descriptions are usually extracted from form factors 
in the limit of vanishing squared four-momentum transfer 
of the probing exchange particle. For example the mean 
square charge radius is defined as the slope of the electric 
form factor at = and the magnetic moment is the 
value of the magnetic form factor at the photon point. 
Although such an approach is suitable to produce pure 
numbers it hardly leads to any insight into the underlying 
structure of the observable. It is for example well known 
how a fermion produces a magnetic moment through both 
its spin and its angular motion, but how does that trans- 
late into the magnetic moment of a bound state e. g. a 
baryon composed of three quarks? 

On the other hand static properties in non-relativistic 
quantum mechanics can be formulated by means of expec- 
tation values involving essentially scalar products of wave 
functions. The nonrelativistic charge radius of a composite 
system of N particles for example is given by 



(1) 



where Qi is the charge of particle i, Xi its position, R the 
center of mass coordinate and Q the net charge of the 
system. A direct relativistic generalization of this expres- 
sion is unknown. In this paper we will focus our attention 



to a quark model description of baryons. The generaliza- 
tion to other systems is then quite obvious. We will show 
that a synthesis of both approaches mentioned before is 
indeed possible - at least if certain restrictions are made 
to the kind of interactions between the constituents of the 
bound system - and leads to new insights into the struc- 
ture of static properties. Moreover the actual computation 
of static moments is then easier and also numerically more 
reliable in comparison to the computation of form factors 
in the hmit — > 0. 

We work within the framework of the Bethe-Salpeter 
equation which has been successfully applied to e. g. baryon 
mass spectra pOElin] and form factors ^E). 

First we briefly outline the Bethe-Salpeter formalism. 
Details may be found in ref. [J. The construction of cur- 
rent matrix elements from Bethe-Salpeter amplitudes is 
addressed, of which a detailed discussion is given in 0]. 
Since by far the most interesting static observables for 
physical applications are mean square charge radii and 
magnetic moments, we show how they can be formulated 
as expectation values with respect to Salpeter amplitudes. 
These amplitudes turn out to be the natural quantities 
which replace the non-relativistic wavefunctions and pos- 
sess a canonical scalar product. Static observables are then 
represented by certain well defined operators whose ex- 
pectation values are computed with the help of this scalar 
product. The emerging operators can be given a natural 
physical interpretation and show very interesting struc- 
tures. To demonstrate the relevance of this formalism, we 
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apply it to a concrete physical model for baryons described 
in refs. [Tl | ^] | ^ to compute nucleon charge radii and mag- 
netic moments and indicate how higher moments can in 
principle be computed. 



2 Bethe-Salpeter equation and current matrix 
elements 

2.1 Bethe-Salpeter equation 

The basic quantity which describes three fermions as a 
bound state is the Bethe-Salpeter amplitude which is de- 
fined in position space through: 

XPa.a.a, (0 ITV^^, (xi (X2 ) 3 (^^s) I P) , (2) 

where ipa i^i) ^'I'e fermion field operators given in the 
Heisenberg picture and Ui arc multi-indices in Dirac-space 
and any internal space which represents a degree of free- 
dom the particle may have. T is the time ordering opera- 
tor. Here |0) denotes the physical i. e. interacting vacuum 
and \P) denotes a three- fermion bound state with total 
four- momentum P on the mass shell i. e. = M^. Be- 
cause of translational invariance it is convenient to intro- 
duce a center-of-mass coordinate X and so-called Jacobi 
coordinates ^ and rj: 



Then by introducing the three-particle propagator 

X Sl,^,,{\P^P^ + \Pr;)S%^,^^,{\P~p,), (7) 

the Bethe-Salpeter equation can be written in a compact 
notation 



(8) 



where a summation over multi-indices and momentum in- 
tegrations is tacitly understood. 

The Bethe-Salpeter equation would be incomplete with- 
out a prescription of how to normalize its solutions. Such 
a prescription can indeed be found (see ref. ^). In a co- 
variant form and using our compact notation it reads; 



p=p 



iA'^^^ + i/\ Xp = 2M^ 

(9) 

where Gq ^ is the inverse three-particle propagator, i. e. 
the inverse of {T)). Here we introduced the adjoint Bethe- 
Salpeter amplitude defined as: 



{P\T{l^lMWaM)i^l.i^^M- (10) 



X := i(xi +2:2-1- 2:3) 

^■.^ Xi- X2 

11 := \{xi +X2~ 22:3) 



X2^X-\i 
X3^X 



2 - 



3'7 



(3) 



The corresponding conjugate momenta are then given by 
the total four-momentum P and the two relative momenta 
pj and Pri : 



P ■■= Pi + P2 + P3 
Pi ^(Pi -P2) 
Pn Upi +P2- 



2P3) 



Pi 

P2 
P3 



if 

3 

3 
3 



Pv 



\Pv 
\Pv 



(4) 



Using this set of coordinates the total momentum depen- 
dence factorizes and we may define the Bethe-Salpeter am- 
plitude in momentum space depending only on the relative 
momenta: 



(27r)4 (27r)4 

xe-'<?'..«)e-'<'^-">Xp.,.,a3fa.P.)- (5) 

One may write the sum of all two-body interactions in the 
form of a three-body interaction kernel: 



if^^^ , , , :=-ft:lo' , .(iP + Pn,Pi.p'f) 

: -1 



-|- cycl. perm.. (6) 



2.2 Salpeter equation 

In order to solve the Bethe-Salpeter equation in physical 
cases relevant for e. g. the structure of hadrons, one ap- 
plies two approximations. First the full propagators are 
replaced by the free ones: 



SpiPi 



(11) 



This approximation accounts for self-energy contributions 
merely by the introduction of effective fermion masses rrii . 
Neglecting retardation effects in the interaction kernels 
leads to the second, so called instantaneous approxima- 
tion. This assumes that there is no dependence of the in- 
teraction kernels on the relative energies in the rest frame 
of the composite system: 



K^^\P,P^,P^,p'^,p',,) 



K^'HlP+Pr„Pi,p'^) 



P={ALO) 



P=(M,0) 



(12) 



W2)(p^,p^) (13) 



(14) 



These conditions can be formulated in any reference frame, 
if all momenta are replaced by: 



P_L :=P- 



p2 



P. 



(15) 
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This space-like vector is perpendicular to the total four- 
momentum and in the rest frame of the system has the 
desired form p± = (0,p). Thus formal covariance of the 
Bcthe Salpeter equation is maintained. 

Adopting both approximations, it is possible to in- 
tegrate out the dependence on the relative energies in 
the Bcthe-Salpetcr equation 0, thus reducing the eight- 
dimensional integral equation to the six-dimensional 
Salpeter equation. This procedure is straightforward if 
there are no two-body interactions in the system. The 
unconnected part of K^^^ however makes the reduction 
more involved. One is then forced to introduce an effec- 
tive three-body kernel which accounts for the effect 
of the two-body interaction approximately (see |21 for de- 
tails) . The effective interaction kernel is then expanded in 

powers of K}^' + , where is the contribution to 



the three-body interaction V^^^ = vj^'^ + V^'^ that couples 

f3) 

to mixed energy components exclusively. Vy[ then corre- 
spondingly is the contribution to the three-particle kernel 
that involves only pure energy components. Up to low- 
est order Born approximation Vj/*"^' the corresponding 
Salpeter equation then reads: 



^MiP€,Pr,) 



A+++ 



A- 



M - n + ie M + Q -ie 
d% d 



7° (g) 7° (g) 7° 



(27r)3 (27r) 

A+++ 



A— 



7° ® 7° 



I 



M - Q + ie AI + n-ie 
terms with cycl. perm, of two-body force. (16) 



Here we introduced the short hand notation A'^^^ 
Af{pi) (g) Af{p2) (g ^13(^3), where Af(pi) are projectors 
onto positive or negative energy respectively and fl := 
wi(pi) -I- 072(^2) + ^3(^3) is the sum of the relativistic 
one-particle energies uji{pi) = + mf. The Salpeter 

equation involves the Salpeter amplitudes, which are pro- 
jected onto purely positive and negative energy compo- 
nents respectively: 



*M(P€-Pr,) := [A+++{p^,Pr,)+A- 



2tt 2tt 



"(P€>Pr,)] 

XM{Pi,P,i) (17) 



The full Bethe-Salpeter amplitude, which is needed to cal- 
culate current matrix elements, can be reconstructed from 
the Salpeter amplitude in the following way: 



XM 



VI 



Go 



where the so called vertex functions F^j were introduced 
which in lowest order in V^f are connected to the Salpeter 



amplitudes by: 

A r A^ 
< = i 



A- 



M -Q M+Q 



7 g) 7 



,0 T^A 
M- 



.7V, 



(19) 



From the normalization condition Q a corresponding nor- 
malization for the Salpeter amplitudes can be deduced, 
which in Born approximation reads: 



{ni\ni) = 



d^p^ d^p^ 
(27r)3 (27r)3 



= 2M. (20) 



Summation over discrete indices is implicitly understood 
here. This norm immediately induces a positive definite 
scalar product: 



(<^l|<?2 



(27r)3 (27r)3 



^*(p|,p^)<?>2(p«,p„), (21) 



whose existence is of utmost importance since static ob- 
servables will be formulated as expectation values with 
respect to this scalar product as announced in the intro- 
duction. 



2.3 Current matrix elements 

To compute any electromagnetic observable, we need to 
know the electromagnetic current (P, A|j'^(x)|P', A') be- 
tween states with total four-momenta P' and P and he- 
licities A' and A respectively, where i^{x) is the current 
operator: 



f{x) =: V(a;)g7^V'(a;) 



(22) 



with the charge operator q. This current matrix element 
can be derived by studying the response of the system in 
an external electromagnetic field in first order of the elec- 
tromagnetic coupling strength jj. One then finds, that the 
corresponding correlation function Gpp, (p^,p^,p^,pjj) 
separates at the poles in the total energy of the system in 
the following way: 



G 



Xp{Pi,Pv) /p|^^.M(o)|p^x XP'{p^,Pv) 



Aujpuj'p, P^ — Up + le 



P'" - LU 



P' + le 



-f- regular terms for P° ujp and P' uip'. (23) 

The Mandelstam formalism and minimal coupling deliver 
an independent way to determine Gpp,. By comparison 
one then finds the following current matrix element 0]: 



{p\rmp') 



(27r)4 J (27r)4 



X rp(Pi,Pv)SF{Pi + iPv) ® Sl{~p(^ + ip,,) 

® S% (P' - p^)rqSl. (P - p,,)P^, {p^ ,pr,), (24) 
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where the adjoint vertex function r]^f{p^,pri) is related to 
^MiP^^Pn) through: 



r^i\pi.Pvh°®l'®l°- (25) 



Note that in the current matrix element H24|) the photon 
couples to the third fermion exclusively. The couplings to 
the other fermions have been accounted for by the factor 
of 3. This is possible, since the vertex functions, which 
describe a composite fermion system, are totally antisym- 
metric. With the explicit boost prescription of the vertex 
function: 

r0{pi,Pv) = SA^^SA^®SA^r^iiA-^Pi,A-^Pr,), (26) 

the time component of the current matrix element in the 
Brcit frame takes the form: 

(ppi/(o)iP)^-3/|i/0n',(p„p,) 

X [Sl,{p^ + ^pr,) ® Sj-i-p^ + p^) ® Sl{M - p,,)] 
X [Sap ®S\^®1\ W®^® fqS%{M - Ap^%„) 

xr^jiA-^\^,A-^%n), (27) 

where V is the space inversion operator, i. e. V{x^,x) = 
(iP , —x). The spatial components arc accordingly: 



{VP\f{Q)\P) 



d'^Pi f d-^Pv -n^ 



(27r)4 J (27r)4^ 
X [SUp^ + ip,,) ® SU-p^+p^) ® S%{M-p^)] 



r^{A-^''p^,A-^\r,). (28) 



Note that in both cases we commuted a triple tensor prod- 
uct of Sap past the three fermion propagators and then 
did an integral transformation to obtain two successive 
boosts. This also explains the appearance of the commu- 
tator [7*,S'yip]„ in the spatial components of the current 
matrix element. 



3 The charge radius 

3.1 From charge distributions to charge radii 

In the next two section we review the definition and precise 
computation of a particular observable namely the charge 
radius of a composite system in the framework of quan- 
tum field theory. The results are well known but we add 
them here for a better understanding. Given some charge 
distribution p{x) one defines its mean square radius by: 



^Jd'x\x\'p{x). 



(29) 



The radius is normalized by the net charge Q, which is 
simply the integral of p{x) over the whole space: 



Q = jd\p{x). 



(30) 



However if the charge distribution has no net charge, the 
normalization 1/Q is of course dropped. If we turn to 
quantum mechanical systems, the charge distribution is 
given by the time component j'^{x) of the four- vector cur- 
rent of the state that describes the system: 



p{x) 



(31) 



Such a state \'ip) can be represented as a superposition of 
momentum eigenstates 



d3p 
UJp 



V'(P)IP). 



(32) 



ip{P) is the wavef unction in momentum space and the 
states \P) are normalized according to: 

{P\P') = 2wp(27r)3,5(3)(P - P'). (33) 

This immediately fixes the normalization of the states \ip)' 



(^1^) = 2(27r)3 / ^^*(P)^{P). 



(34) 



Let us further investigate the charge distribution by in- 
serting ((221) into (pTT|l : 

1 f (PP f d^P' , 
p{x) = ^— / / exp (i(P - P') • a;) 

ivm J J ^p' 

xripmp'mfmp'). (35) 

We used space translation invariance here to separate the 
spatial dependence. As is well known the integral / d%: exp(ip- 
a;) is a representation of the delta distribution. So it fol- 
lows 



J d='x|a;|2exp(i {P - P') ■ x) 



(2nf 



{Vp-Vp')H'^^\P-P'). (36) 



Using 1123 and we obtain: 

4 J UJp J UJp' 

X exp (i {up ~ ujp')x°) i,*{P)ij{P') 

X (P|j"(0)|P')(Vp- VpO'<5'^'(P-P')- (37) 

On the right hand side of this equation we may now inte- 
grate by parts twice and subsequently do one of the two 
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momentum integrations: 

i^(p)p 



[(Vp-Vp,f(P|j"(0)|P')]p'=p 



(Vp - Vp' 



LUp LUp' 



(P|/(0)|P) 



(Vp - Vp,) 



LUp UJpi 



P'=P 



[(Vp-VpO(P|/(0)|P')]^, 



p'=p 



(38) 



The last of this three terms vanishes, because Vp and 
Vp' change sign under space reflection, in other words 
are of odd parity. So if we assume that the states \P) have 
definite parity then (P|j°(0) Vp|P) = 0. 

So far we have considered wave packets that consist 
of a superposition of states with different momenta. To 
obtain states with definite i. e. sharp momenta consider 
first a Gaussian wave packet with a width proportional to 
some parameter A: 



V<P)=exp(-|P|V(2A2)) 



(39) 



Before we let A go to zero to define a definite momentum 
state let us inspect the second term in (|38ll . Because for 
the Gaussian wave packet from (|39|l the wavefunction is 
real i. e. ip*{P) = ip{P) we have: 



(Vp- VpO 



LUp 



Jp, 



0. 



p'=p 



(40) 



Therefore also this term does not contribute to the charge 
radius and we are left with the first term in H38(l only. 

Let us now turn to the limit A — > again. Since in this 
limit exp(— |Pp/A^) is another representation of the delta 
distribution we find 



Lop{i,{pmp)) 



J(3)(p) 



(41) 



2(27r)3 ■ 

Inserting this together with the first term of H38|l into the 
basic definition of the charge radius (^5)) and performing 
the final momentum integration yields 



8A/g 



(Vp-VpO^(Pb-"(o)|p') 



P'=P=0 



(42) 



where M is the rest mass of the system. The current ma- 
trix element appearing here is given in the Breit frame if 
we make the following transformation: 



(Vp - VpO^(P|/(0)|P') 



P'=P=0 



Ap{VP\3^Q)\P) 



p=o 



We then finally end up with the expression: 
1 



8MQ 



Ap{rp\fmp) 



(43) 



(44) 



p=o 



3.2 From form factors to charge radii 

In the last section the derivation of the charge radius op- 
erator started from defining the mean square radius of a 
charge distribution. As is well known there is another def- 
inition of the charge radius that involves the electric form 
factor of the system. There the charge radius is defined as 
the slope of the electric form factor at the photon point. In 
this subsection we investigate the interconnection between 
both definitions and show that they indeed coincide. 

Let us briefly recall some basic definitions in the con- 
text of form factors. From current conservation and Lorentz 
invariance the electromagnetic vector current of a spin- 1/2 
state can be parameterized as follows: 



(P',A'|/'(0)|P,A) =e^2v(^') 



2M 



ux{P). 
(45) 



Pi and P2 are the Dirac and Pauli form factors respec- 
tively. The Dirac form factor is normalized to the charge Q 
whereas the Pauli form factor is normalized to the anoma- 
lous magnetic moment k. Both form factors are functions 
of the squared invariant momentum transfer := —q^ = 
— (P' — P)^. The Dirac spinors are normalized in a Lorentz 
invariant fashion: 

ux'{P)ux{P)^2MSyx. (46) 
Using this normalization one shows that 

uy {VP)ux{P) = 2 + 0V4Ja'a (47) 

as well as 

iiA'(PP)7°«A(P) = 2A/<5a'a. (48) 

Using both expressions one can write the time component 
of the electromagnetic vector current (|45|) in the Breit 
frame as 



(PP,A|i°(0)|P,A) =2eMGi5(g2 



(49) 



where Ge{Q^) is the electric Sachs form factor. It is de- 
fined together with the magnetic Sachs form factor as a 
combination of the Dirac and Pauli form factors: 



Gb(Q2) :=Pi(g2)_ 
Gm(Q') :=Pi(g') + 



4A/2 



F2iQ' 



(50) 
(51) 



The mean square charge radius is defined as the slope of 
the electric form factor at the photon point i. e. at = 0: 



r- = -- 



6 dGB(g2) 



Gb(o) dg2 



(52) 



Since G£;(0) = Pi(0) = Q the normalization 1/G£;(0) is 
dropped in case the net charge vanishes. From this defini- 
tion together with H49() we then find 



_3 d_ 

Mgdg^ 



(PP,A|/(0)|P,A) 



(53) 



^=0 
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This result has to be compared to the one that we obtained 
in the previous section, namely equation (|44() . There we 
found the Laplace operator with respect to P instead of 
a single derivative with respect to acting on the cur- 
rent matrix element. However both expressions turn out 
to be exactly equal: From the parameterization of the vec- 
tor current l|45() it is clear that the current matrix element 
("PP, A| j" (0) I P, A) depends on Q"^ . In the Breit frame the 
dependence of on the momenta of the incoming and 
outgoing bound states becomes rather simple. It reads 
(52 _ 4|p|2^ Then for any function / depending on 4|Pp 
the following identity holds: 

Z\p/(4|PH = 4 {Ap\P\^) 4rif^Q^) = 24;^/(Q^). 



dQ2 



dQ2 



(54) 

Inserting this into (|44|l we see that it coincides with (|53|) . 
Thus the definition of the charge radius from form factors 
is exactly equivalent to that from charge distributions. It 
must be noted however that expression H44|l is somewhat 
more general than (|53|) in the sense that it is valid for 
particles with arbitrary spin. Nevertheless one can show 
that both expressions lead to the same result. We decided 
however to start from H53|) simply because it contains only 
a first order derivative. 



3.3 The charge radius as an expectation value with 
respect to Salpeter amplitudes 

We now want to analyze the charge radius in the Bethe- 
Salpeter framework for a three-quark system. It is essen- 
tial to know the Q^-dependence of the current matrix el- 
ement (|27|l . So let us inspect its P-dependent part alone: 



C.2 



>SUM-Ap'%r,) 



S'a, 



l]fiAp''p^,Ap' 



>,). (55) 



We now exploit an important property of Lie groups, 
namely that every group element may be represented as 
an exponential mapping of the Lie algebra: 



l\f{A-p'%^,Ap''p,,) 

= exp(-2i77(P).K)/(p5,p„) (56) 



The parameter rj, commonly called rapidity, is defined as 
follows: 

where the last equality follows from Breit frame kinemat- 
ics. The operator K is an infinitesimal boost. The genera- 
tors of the Lorentz group are given by the following skew 
symmetric tensors (see e. g. [HI): 



(58) 
(59) 



Because of skewness there are six independent quantities. 

are the three generators of boosts and the remaining 
three operators generate rotations (in fact they are the an- 
gular momentum operators). In momentum space we have 
Jp" = i{p^d/dpi, — p^d/dpfj). 5*^" are the corresponding 

generators in Dirac space. The infinitesimal boost K then 
simply reads 



jOi 



Pi 



'Oi 



1 



Pew, 



Pa 



„0 d 



PridK, P'tdp^ 



(60) 



Inserting H56(l back into the current matrix element (|27|l 
we find 



(PP|/(0)|P) = -3 



Pi 



(27r)4 J (27r)4 
X exp(-2ir7(P) -K) [1(g) 1(g) -f°qSl{M ~ p^)] 

xr^ip^.p,). (61) 

Since the charge radius is proportional to the slope of the 
current matrix element (|61|l at = we are interested in 
the term of the expansion linear in and thus - because 
|r7(P)p is of order - Hnear in \t]{P)\'^. Writing the 
expansion explicitly out up to this order we have: 



(PP|j°(0)|P) 



-3 



2MQ 



d^Pr, 



(27r)4 J (27r)4 
X [S\.{jp^ + Ip^)® Sl{~p^ +P^)g) Sl{M - p^)] 



3 

E 



if{P)ri\P)k^k^ 



X [l®I®7VI(M-p,)]P/,(p^,p,,) + 0(r74). (62) 

By inserting this into (|53|l . the charge radius then takes 
the form: 



18 



3 



Q2=0 



a Pi 
(27r)4 



Pn -r^ I \ 
l^MyPi^Pry) 



(2^)4' 

X \s\{p^ + \pr,) ® sU-p^+Pr,) <E> sUm-p,,)] k'k^ 

x[lg>lg> j°qSUM - p,,)] P/,(pc,p„). (63) 



The integration over the relative energies can now be per- 
formed by using the partial fractions decomposition of the 
fermion propagators: 



sUp.) 



AtiPi) 



A- (Pi) 



p°-Wi(pi)+ie p'^ + LJi{pi) - ie ^ 



(64) 
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With the aid of Cauchys theorem both integrations can 
be performed and one obtains: 



3 



18 



d 
dQ2 



(27r)3 J (27r)3 

X k'^k'^q^ 



A+++ 
{M - i7) " 
A— 



A- 



(M + Q) 



(M - [2) 

X [7" «) 7° 



(A/ + i7) _ 
'l'']r^APi,Pv)- (65) 



Now (j3 denotes the charge operator acting on the third 
fermion. After integration the boost becomes: 



K' 



2W2 ' 



(66) 



Note that we used the anticommutator {7'^, a*}+ = and 
the relation Afa^ = Af ±p{/uJi here. Now by using re- 
lations (|19|l and H25|l one can replace the vertex functions 
in l|65|l by Salpeter amplitudes. The result is an expec- 
tation value with respect to the Salpeter scalar product 



18 



3 

MQ 



d 

dQ 



X {i>Uk'^k'^f\<i>i,) 



(67) 



Since K' is a tensor operator of rank 1, that is a vector 
operator, k'^k'^ is a Cartesian tensor operator of rank 
2. As is well known, every Cartesian tensor may be de- 
composed into irreducible representations of the rotation 
group 5*0(3). The decomposition of a rank 2 tensor is 
given by: 

T,, = \^{T)5,, + \{T,,-T,i) 

+ i(r,j+r,,-|tr(T)<5,,). (68) 

According to their transformation properties under rota- 
tions, the first term belongs to the scalar representation, 
the second to the vector representation and the last to 
the five dimensional representation of spin 2. Let us now 
address the question, which of these representations will 
vanish due to selection rules in the scalar product in (|67|) . 
Let us start with the vector representation: 



3 

E 



d 
dQ 



{<PU\(k'^k'= -k'Ui'^)(t\i>i,) 



0. (69) 



This is so, because rj'' {P)rj^ [P) is symmetric, whereas 
k''^k'^ — k'^ k'^ is antisymmetric under the exchange of 
indices. For the spin 2 representation we cite the Wigner- 
Eckart theorem and in particular the triangularity relation 



which states, that for spherical tensor operators of rank 
k: 



{J 



1 iTf'l b2> = unless | ji - 32] < k < n + .h- (70) 



In our case ji = j2 ~ \ and A: = 2, so the spin 2 rep- 
resentation in (|68|1 gives no contribution. Only the scalar 
representation contributes and we get from and l|()7jl : 



6 



MQ 



;v'iP) 



Recalling the definition of the rapidity (|57|l we find: 



dQ- 



rV\P) 



4M2' 



which brings us almost to our final result 
3 



2A/Q 



M2 



(71) 



(72) 



(73) 



By rewriting k' in terms of one-particle coordinates we 
find: 



K" = 



1 



n 



h. C. 



It is also useful to define: 



1 ^ 

-R := ^ E '^"^^Pc 

ct=l 



(74) 



(75) 



The expression (|73|) is still not symmetric in all three par- 
ticles. The third fermion seems to play a special role. How- 
ever this asymmetry is only due to the fact that in deriving 
the current matrix element we exploited the total asym- 
metry of the vertex functions under particle interchange 
and coupled the photon to the third fermion exclusively 
accounting for the other couplings by multiplying with a 
factor of 3 (see 0). If we reverse this procedure, cancel the 
factor of 3 and symmetrize the expression over the three 
particles, we end up with a symmetric version: 



1 



Q = l 



n 
Ji 



R 



h. c. 
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(76) 



Now the sum runs over all three fermions. Note that we 
also used the norm of the Salpeter amplitudes (|20|l to re- 
place 2M by (^-^.K). 



3.4 Interpretation 

Having derived an analytic expression for the mean square 
charge radius of a bound three-fermion system with in- 
stantaneous interaction kernels, it is worthwhile to give 
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fcrinioii 1 



*1 



\_ccntcr of mass 



' fori 



4 The magnetic moment 

4.1 The magnetic moment as an expectation value 

To derive magnetic nioments from form factors in an anal- 
ogous way we start from the parameterization (|45|) of the 
electromagnetic vector current: 

(7'P,A'|j + (0)|P,A) 

= e [^^i(g') +^^2(Q')] u),,{rPh+ux{P), (81) 

'fcmiioii 3 where the "+" -component of the current is defined by: 

j + (0)=/(0)+i/(0). (82) 

Note that this definition deviates from the definition of 
the components of a spherical tensor operator of rank 1 
which are usually given by: 



Fig. 1. O denotes the origin of the reference frame, R the 
position vector of the relativistic center of mass and xi through 
X3 the position vectors of the three fermions. 



the result a meaningful physical interpretation. To inter- 
pret the operator between the Salpeter amplitudes it is 
useful to note that iVp^ is the position operator in mo- 
mentum space: 



— 



(77) 



Consequently the quantity R as defined in (|75(l is the 
canonical relativistic center of mass of a three-particle sys- 
tem: 

1 ^ 

^=-^tj„*„. (78) 

a=l 

At fermion momenta small compared to their masses, i. e. 
\Pa\ ^ ™ai we have uja "Tia and thus the expression 
reduces to the well known nonrclativistic center of mass: 



-Rnr 

The expression 



1 



mi + TO2 + ^3 



ITlai 



(79) 



a=l 



1 ^ 



(80) 



13=1 



then corresponds to the position of particle a as measured 
from the relativistic center of mass. Fig. ^ illustrates the 
situation. Since iVp^ — R in the difference between two 
vector operators, it is invariant under translations and 
consequently the mean square radius H76|l is translation- 
ally invariant. Finally we want to call attention to the 
relativistic factor fl/M in H76|l which weights the relative 
distance of each fermion with the collective relativistic en- 
ergy. Therefore an enhancement can be seen the more rel- 
ativistic the system is. 



T. 



[1] 



T^(ri±iT2) and T^'' :- T3. (83) 



The total spin of the system makes a spin fiip of course so 
we have A' 7^ A. The spin polarizations will be fixed later. 
Evaluation of the spinorial part of this equation yields: 



Sa'(7'P)7+"a(P) -2^0^ 



(84) 



Together with the definition of the magnetic Sachs form 
factor H5U|I we then get from H81|) the relation: 



Gm(Q') 



(7^P,y|j+(0)|P,A) 



(85) 



which expresses the magnetic form factor in terms of spa- 
tial components of the current matrix clement. The mag- 
netic moment is defined as the value of the magnetic form 
factor at the photon point: 



(m> :=Gm(Q' = 0). 



(86) 



Because of the denominator in (|85|l taking this limit re- 
quires some care. We need to know the Q^-dependence 
of the current matrix element. But first let us choose the 
three-momentum transfer to point in the 3-direction from 
now on: 

g:= j = v/Q^ea. (87) 

With this choice we have [y^ ^ Sap\- = such that the 
current matrix element H28(l now becomes: 



{vp\3+mp) 



(27r)4 J (27r)4^ 

X [Slip^ + \Pr,)® SU-p^+p^) ® SUM-Pr,)] 



'St 



I (g) I (g) qj+Sj,{M - Ap' p^) 
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To extract the Q^-dependence of this expression we ex- 
press the boost as an exponential just hke in the previous 
section on charge radii: 



(7'P,A'|j+(0)|P,A) 



(27r)4 J (27r)4 



X ^M,A'(P€,Pr7) 

X exp(-2ij7(P) • k) [I® I® g7+S'|,(Af -p,,)] 

X P//,a(pc,P,)- (89) 

This time however the boost generator also acts on the 
Dirac space of the third fermion as can be seen from the 
current matrix element (|88|l : 



„0 _d i _d _ _d i _d 

^(.dpl ^^Wl Pvdplj Pnd^ 

-I- i (a* (8) I (8) I + cycl. perm. 



(90) 



Inserting the expansion (|89|l into H85(l and taking the limit 
^ then shows that terms with 0{ri) > 1 vanish 
because ri{P) is of order \/Q^- Concerning the first or- 
der term we find with the special choice of the three- 
momentum transfer H87() and the definition of the rapidity 



lim 1 = lim \Q __ — ^ 



63 



2M 



(91) 

Therefore the static limit can safely be taken and we find 
for the magnetic moment: 



d'^P£, f d'^prj —A 



4A/2 J (27r)4 J (27r)4" 

X ii^3 [I ® I ® qi+sUM - Pr,)] r^iAPe^Pn). (92) 

We inserted a factor 1/2 A/ in this expression since the 
wavefunctions are normalized to 2Af as can be seen from 
(|20() . Integration over the relative energies can now be 
done after replacing the fermion propagators by their par- 
tial fraction decomposition Ht)4|) : 



(a*) = — - 



d^Pi f d^p^ -A 



(27r)3 



(27r)3' 
A- 



{M - n) {M + Q) 

A- — 



{M - i7) {M + Q) 



[7"®7°®7°]A1,a(P€,P^ 



(93) 



with the tensor operator: 



-h. c. 



1 



7 

Pi 



2M [ 2^3 
f2 



\ (cji - tJa) ig^ + (^1 + ^2) ijr 



2t^3 



1(g) II(g)ia*a-' 



2^1 



^2 I , 



(94) 



where we also used the anticommutator {7^,0;*}+ = 
and the relation Ai^a? — Af ± p^/tUi. Note that the 
-component in the second index of P*-* in H93|l has 
to be taken in the sense of l|82(l . Before we analyze this 
expression further let us replace the vertex functions in 
by using the relations l(T^ and ipSjl to arrive at the 
compact notation: 



2M 



(95) 



Since P'-* is a product of two vector operators it consti- 
tutes a Cartesian tensor operator of rank 2, which can be 
decomposed into irreducible representations of the rota- 
tion group according to H68|) . Just as we did when deriving 
the charge radius, we may show that the contribution of 
certain representations vanish. The scalar representation 
gives no contribution because of the m-selection rule of 
the Wigner-Eckart theorem, which states that 

(ji,TOi|Pj''l|j2,nT-2) = unless mi - m2 = q. (96) 

In our case mi = i, 7712 = ~^ and q = 0. The spin 2 rep- 
resentation vanishes because of the triangularity relation: 



'1 i|r[2i|i 

\ 2 ' 2 I g I 2 ' 



0. 



(97) 



We are thus left with the antisymmetric representation be- 
longing to spin 1 which we may write as a vector product: 



p3+m ^ (^31 ^^^32^1] ^^J2 ^+,kF''' 

3,k=l 



(98) 



where it is understood to take the spherical "-|-l"-component 
of the vector product as defined in H83|) . Note that since 
F"^^ is contracted with the skew tensor tijk, the last term 
in H94|) that is proportional to p|P3 vanishes. Inserting (|98|l 
back into (|95() and choosing the spin projections ^' = \ 
and A = — i then yields: 



(m> = ^«i/2l^ E ^+3kF"'^^\'^ 



A/, -1/2/ 



(99) 



By using the Wigner-Eckart theorem once again we re- 
move the spin-fiip and turn the expression in an expecta- 
tion value: 



3 
2M 



«i/2l E ^3,fcP^'=g3Ki/2), (100) 
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To simplify H94(l further we replace the relative coor- 
dinates by one-particle coordinates: 




h. C 



1 ^ 

\_d i \^ , , ; 8 



2^3 



n 

a=l 



(101) 



Since in the expectation value (|100|l F^*"' is contracted with 
the skew symmetric tensor eijk it is suggested to define: 



(102) 



^Ra obviously the total angular momentum of the three- 
quark system with the correct center of mass motion re- 
moved. Furthermore we identify the spin operator S = 
il7 in the following contraction: 



3 



2i 



We then have: 



<7^ I 



2iS\ 



(103) 



3 



j k 



f 3 



h. c 



(104) 

This expression is still not symmetric in the three fermions, 
so in the final step we symmetrize over the three fermions 
in the same way as we did already when deriving the 
charge radius: 



Km 



Ml 



(105) 



where we defined the magnetic moment operator (1 which 
follows from symmetrizing H1U4|) : 



, 1 



(106) 



with the one-particle spin operators: 

:= I!/2(E)l(g)l 

52 ■■= I® i:/2(8) I 

53 := l(g)l(E)S/2. 



(107) 



4.2 Interpretation 

As has already been shown in the interpretation of the 
charge radius, the term H8U|1 corresponds to the position of 
particle a as measured from the center of mass of the sys- 
tem. One is thus naturally led to interpret defined in 
1)102(1 as the angular momentum (operator) observed from 



the relativistic center of mass. As already mentioned Si, 
S2 and ^3 are one-particle spin operators. We therefore 
conclude that the magnetic moment of the system can be 
decomposed in contributions of the fermion angular mo- 
menta and their spins: 



(108) 



with (fJ-L) being the contribution of the angular momenta 
of the three fermions: 



(mi 



1 



9 I M ^ 



3 

Ra 



and {^s) the contribution of the fermion spins: 



(a^s) 



(109) 



(110) 



Such a decomposition into spin and angular momentum 
contributions is not possible when extracting the magnetic 
moment from a form factor. It is thus another benefit of 
the approach to static properties presented in this work. In 
(|106() we discover the same relativistic weight factor Q/M 
as has already been found in the charge radius. When tak- 
ing the nonrelativistic limit, the operator jl (|106|) becomes: 



2mn 



^H^Sl (111) 



a=l 



2ma 



Except for the center of mass correction this expression 
is well known. We are thus led to conclude that we have 
found the relativistic generalization of the nonrelativistic 
magnetic moment operator. 

Although the expression H86|l we started from holds 
only for spin 1/2 states, one can, as in the case of the 
charge radius, generalize the formalism to arbitrary spins 
and the result is not different from the spin 1/2 case. 



5 Extension to higher moments 

The formalism presented so far paves the way for the cal- 
culation of higher moments as well - an issue that we 
briefly want to touch upon in this section. We take the 
electric form factor as an example and work accordingly 
with the "time" -component of the current matrix element 
(|27|l . An arbitrary moment (m) of a charge distribution is 
then given in general by: 



E 

il ,S2 , ■■■1*71 = 1 



Oil i 



XXX 



(112) 
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where Oi 



is a tensor of rank n. which depends on 



the moment to be computed. For example for the charge 
radius, considered so far, O is simply: 



1 . 



(113) 



By similar steps leading from eq. to eq. (|^ we get: 



2M 



E 



_d d_ 



d 



dP 



-(^P|/(0)|P) 



(114) 



The current matrix element appearing here is defined in 
eq. H27() . As before its P-dependent part is given by an 
exponential of infinitesimal boosts as in eq. (|56|l . Because 



lim riiP) = and ^—-nHP) 



f =0 



(115) 



we find: 
d d 



d 



Qpn Qpi2 Qpi 



■ exp(-2ir/(P) • K) 



t3L 



p=o 



(116) 



Note that to every from our starting equation (|112|l 
corresponds now a boost generator K^. Using this result 
we get from eq. pi4f) : 



1 ^ 

9 M ^ 



2M 

x(-3) 



2l ,l2,...,in — 1 



d\—A 



A,f(P€,P»7) 



(27r)4 7 (27r)4 

X [5].(p5 + ip„) (g) S'|^(-Pc + p^) «) S%{M - p,,)] 
1 



A/" 



X []I®II®7Vf(A^-P^)] A^/(PC,P^)- (117) 

Integrating out the dependence on the relative energies 
after replacing the propagators according to eq. H64|) then 
results in: 



E 



1 



X > On......„ (Klj^KK ■ ■ ■ Kiun 

+ off-diagonal matrix elements (118) 



where K'^ is defined in eq. (|66|l . For n > 2 we also find 
terms involving matrix elements between different energy 



components of the vertex function, i. e. between the sub- 
spaces of purely positive and negative energy components 
(denoted "off-diagonal matrix elements" in eq. (|118|l ). Un- 
fortunately these terms cannot be expressed in a generic 
way and have to be calculated explicitly for the moment 
under consideration. One might however expect that these 
additional contributions arc in fact small; first because 
the negative energy components correspond to the "small" 
components of the Dirac equation and thus vanish in the 
non-relativistic limit and second because both energy sub- 
spaces are orthogonal. Note that although the first term 
of eq. H118|l also involves matrix elements between differ- 
ent energy subspaces of the Salpeter amplitudes, one can 
show that these do in fact vanish. 

Finally we may symmetrize the expectation value in 
eq. H118() over the three fermions to obtain: 



1 



E 



o. 



K/lE^^''U^^U---^.'lo9c.|<Z': 



off-diagonal matrix elements, (119) 



where K'/^^ is defined as: 



K" 



n 

M 



R 



(120) 



If e. g. we insert eq. (|113|) into the final result (|119|l we 
instantly obtain the charge radius expression (|76|l . What 
has been said about its interpretation also applies to eq. 
(|119() in its general form. 

In this sense a generalization of the formalism pre- 
sented in this work to arbitrary moments is possible, al- 
though those discussed in detail, namely the charge radius 
and the magnetic moment are by far the most interesting, 
having in addition the soundest empirical basis. 



6 Application to static properties of baryons 

We would like to illustrate the relevance of the preceeding 
theoretical considerations by applying them to an exist- 
ing physical model. In refs. PEIE] a relativistic covariant 
quark model for baryons is treated, based on assumptions 
which also entered the work at hand, i. e. instantaneous in- 
teraction kernels and free fermion propagators correspond- 
ing to effective fermion masses. The model successfully de- 
scribes mass spectra of strange and non-strange baryons 
up to the highest orbital and radial excitations employ- 
ing a linear confinement potential and a residual inter- 
action based on an effective instanton force. The seven 
parameters entering the model are fixed by a fit to the 
best established resonances. We use the wavef unctions, 
i. e. Salpeter amplitudes, that have been obtained by solv- 
ing the Salpeter equation within this model to compute 
the expectation values of the charge radius and magnetic 
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moment operator derived in this work. Since no further 
parameters are introduced the results are genuine predic- 
tions. 



6.1 Nucleon charge radii 

The proton charge radius that we obtain by computation 
of the expectation value (|76(l amounts to 

\/ \ /proton 

0.86 fm, (121) 

in excellent agreement with the experimental value of 
0.87 ± 0.008 fm from ref. |7]. The mean square charge ra- 
dius of the neutron however results in 

(r2)ncutro„- -0.206 fm^ (122) 

and overestimates the empirical number of 
-0.1161 ± 0.0022 fm^ from ref. [3 by 77%. Within the 
same model the authors of ref. jj] have calculated the 
neutron electric form factor and extracted a mean square 
charge radius of — O.llfm^ from it. The procedure how- 
ever was numerically erroneous and a reanalysis, improv- 
ing the numerical precision, resulted in a radius that is 
indeed compatible with our result. 



hyperon 


experiment jj] 


this calculation 




[M/Miv] 


[^/Miv] 


A 


-0.613 ±0.004 


-0.61 


S+ 


2.458 ±0.01 


2.51 


E° 




0.75 


E- 


-1.16 ±0.025 


-1.02 


^0 


-1.25 ±0.014 


-1.33 




-0.6507 ± 0.0025 


-0.56 



Table 1. Hyperon magnetic moments compared to the empir- 
ical values. 



nucleon resonance 


h 


magnetic moment 


P"(1440) 


1/2 


1.55 




-1/2 


-0.98 


S"(1535) 


1/2 


0.37 




-1/2 


-0.1 



Table 2. Prediction of magnetic moments of selected excited 
nucleon states. Iz means third isospin component. 



6.2 Baryon octet magnetic moments 

In the same model we have computed the nucleon mag- 
netic moments using our formula (|105|) . For the proton we 
find a magnetic moment of 

(Ai>proton = 2.77 ^AT (123) 

in perfect agreement with the empirical value of 2.793 /iat. 
The magnitude of the neutron magnetic moment 

(Ai)noutron = -1-71MW (124) 

is rather small if compared to the experimental magnetic 
moment of —1.913 yUAr. 

In addition to the nucleon magnetic moments we have 
also calculated those of the strange octet baryons because 
they are experimentally well covered. Table ^ compares 
our results to the empirical values. The results are in ex- 
cellent agreement with experiment. The largest deviation 
of 14% is seen with the S'^-magnetic moment. 

Since efforts are being made to measure also magnetic 
moments of excited nucleon states like the S^^(1535) as 
mentioned in ref. |H|, we contribute some selected predic- 
tions here. The magnetic moments of the nucleon Roper 
resonance (P^^(1440)) and the lowest lying state with to- 
tal spin 1/2 and negative parity (8^^(1535)) are shown in 
table 13 The formalism allows the computation of mag- 
netic moments of baryons with arbitrary spins and their 
radial excitations which will be the subject of a subsequent 
publication. The same is true of course for the charge ra- 
dius. 



As has already been indicated, the magnetic moment 
may be decomposed in spin and angular momentum con- 
tributions according to eq. 11081 This decomposition en- 
ables us to carry out a numerical analysis of the magni- 
tudes of both spin and angular momentum contributions. 
Note that such a study is not possible by relying on form 
factor calculations because there only the total magnitude 
of the magnetic moment can be extracted. Table|311ists the 
contributions of spin and angular momentum to the mag- 
netic moments of proton and neutron. The analysis shows. 





2{/^s> 
b/^iv] 


2(ms) 
[%] 


(/iL> 


(Mi) 

(m> 

[%] 


proton 


2.53 


91 


0.24 


9 


neutron 


-1.59 


93 


-0.12 


7 



Table 3. Contributions of quark spins {2{fis}) and angular 
momentum to the net magnetic moments of proton and 

neutron. 



that the contribution of the quark spins exceeds the contri- 
bution of the quark angular momenta by far. One can state 
that a good 90% of the magnetic moment is coming from 
quark spins which is due to the fact that the quarks are 
dominantly in a relative S'-wave. This result also explains 
in part the success of the nonrelativistic quark model in 
predicting the magnetic moments. Our analysis shows that 
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50 100 150 200 250 300 350 
quark mass [MeV] 

Fig. 2. Fraction of the total magnetic moment carried by quark 
spin and angular momentum respectively of proton and neu- 
tron respectively as a function of the quark mass. 



by neglecting the angular motion of the quarks by assum- 
ing that the quarks are in a relative 5'-wave, the induced 
error is in the percent region. We should mention that for 
the 8^^(1535) the absolute value of the spin contribution 
is only a quarter of the angular momentum contribution 
and opposite in sign. Since this resonance is dominantly a 
P-wave the spin has to be aligned antiparallcl to the an- 
gular momentum to result in a state with total spin 1/2. 
The preceeding discussion is however only true if we work 
with a constituent quark mass of 330 MeV. 

Wc may however carry this analysis further by study- 
ing the evolution of the spin/angular momentum distri- 
bution with smaller quark masses. Note that the quark 
model described in refs. m^lEl assumes isospin symmetry 
between up- and down-quark and thus there is only one 
mass parameter for the nucleon. At different magnitudes 
of this mass parameter we have now fitted the remaining 
six parameters of the model to the baryon spectra. We 
might of course not expect to reproduce the spectra as well 
as with the original value of 330 MeV but at least we were 
able to keep the ground states i. e. the nucleon and the A- 
particle at the empirical values. We achieved a quark mass 
as small as 25 MeV before numerical restrictions impeded 
us to go any further. Figure 16.21 shows the effect on the 
spin and angular momentum contribution to the magnetic 
moment of proton and neutron. We see an almost linear 
decrease of the spin contribution from its original value of 
a good 90 % at 330 MeV to roughly 60 % at 25 MeV. At 
the same time the angular momentum contribution gains 
in magnitude correspondingly to roughly 40%. Although 
one looses the concept of constituent quarks at such small 
quark masses, the analysis nevertheless shows that when 
approaching the chiral limit, spin and angular momentum 
contribution to the magnetic moment become of the same 
order of magnitude. 



7 Conclusion 

We have shown how the charge radius and the magnetic 
moment of a bound three-fermion system with instanta- 
neous interactions can be formulated as expectation val- 
ues with respect to Salpeter amplitudes. The correspond- 
ing operators turned out to be natural relativistic gener- 
alizations of their non-relativistic counterparts. We also 
indicated how the formalism may be extended to higher 
moments as well. A first application of the formalism to a 
relativistic quark model for baryons with instantaneous 
interactions described in refs. [T]P1|^ results in a good 
description of the nucleon charge radii and baryon octet 
magnetic moments except for the neutron radius. Predic- 
tions have been made for the magnetic moments of the 
P^^(1440) and 8^^(1535) resonances. Wc found in addi- 
tion an interesting dependence of the nucleon magnetic 
moments on the quark masses. In particular, when con- 
stituent quark masses are decreased to values almost as 
small as current masses, the spin contributions become 
equal in magnitude in comparison to the contributions of 
internal angular momenta. Static observables of systems 
with spins other than 1/2 like e. g. the baryon decuplet 
will be studied in the future. 
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